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ABSTRACT 



ANALYSIS OP INELASTIC ALPHA PARTICLE SCATTERING 
USING DIRECT INTERACTION THEORY 



by 



Herbert Prank Swanson Jr, 

The theory of the direct interaction is discussed 
and its applicability to inelastic alpha particle 
scattering to low lying states is shown. Two results 
of this theory are developed; (1) the diffraction model 
of Blair and (2) the distorted wave Born approximation, 
where the phenomenological approximation to the optical 
wave function suggest ed.^by McCarthy and Pursey is used. 
Pits to Ar^°( ©C, o£') Ar data are shown for 32.8 and 
41.0 Mev(Lab) alpha particles. In the diffraction model 
analysis, R Q =* 6.45 Fermis for both energies. 

/^R c — .82 Fermis for the lower energy data and .71 
Fermis for the higher energy data. /^?R 0 =* *67 and .56 
Permis for the lower and higher energy, respectively. 
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CHAPTER X 



THE DIRECT INTERACTION 

There are two basic mechanisms by which a nuclear 
reaction can proceed; (1) the direct interaction, and 
(2) the formation of a compound nucleus. The theoretical 
angular distribution predicted by- each mechanism is 
different. To see the nature of these two processes, 
consider the incident particle to be a nucleon. Thus 
the independent-particle shell model is a good 
approximation to the nuclear structure. Here the many 
body problem is replaced by an effective single particle 
potential which describes the force acting on each 
particle, this potential being different for different 
shell model configurations. 

As a particle enters the nucleus, it experiences 
a potential described by the other A nucleons in their 
ground state. configuration. The projectile may then 
simply enter the nucleus, be deflected by the nuclear 
potential, and emerge again at a different angle but 
with the same energy in the center of mass system. 

This is called direct elastic scattering. If a 
reaction ensues, the energy transferred to the nucleus 
is determined by the interaction of the incoming 
nucleon with one in some shell configuration and by 
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the nuclear density. This later fact modifies the 
problem as neither of the colliding particles can have 
an energy below that of the ground state energy of the 
nucleus. This is because the hole left by the struck 
nucleon is the only level available below the ground 
state energy. The case where the incident particle 
replaces the struck nucleon is called exchange scattering 
and will not be considered here. Those nucleons excited 
to high level states should be described by a different 
potential than those excited to low lying states, 
hence the effective potential is a function of the 
bombarding energy. If either of the colliding nucleons 
has an energy greater than its separation energy, it 
may leaye the nucleus without further reaction save 
deflection by the nuclear potential. This is described 
as a direct interaction and usually occurs in a time 
comparable with the time it takes the projectile to 
traverse the nucleus. If, on the other hand, neither 
of the nuclei has energy sufficient for separation, 
both will undergo further collisions eventially 
spreading their energy over the whole nucleus. For a 
time, the nuclear state will grow increasingly complex 
until it reaches a statistical equilibrium. This state 
will be a complicated admixture of shell-model and 
collective configurations involving excitation of the 
many degrees of freedom of the nucleus. Eventually 
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enough energy will be concentrated in a particle or 
group of particles to allow escape from the nucleus 
and upon emission the nucleus will decay into some 
new excited state. This process is referred to as 
scattering via the formation of a compound nucleus 
and is usually characterized by the long time necessary 
for its completion, which is about a thousand times 
that for the direct interaction. Since these two 
processes exhibit widely different cross sections, a 
method for separation into the direct interaction 
and compound nucleus components is desirable. No 
suitable technique, however, has been found. Instead 
reactions are chosen which are predominantly either 
direct interaction or compound nucleus types and where 
errors made in neglecting the other are small. 

Direct interaction processes involve only a few 
of the many degrees of freedom of the nucleus. The 
minimum number, in fact, would be just those necessary 
to specify the ingoing and outgoing channels of the 
reaction. The concept of the direct interaction can be 
extended to any projectile if it excites only one degree 
of freedom in the nucleus specifically a collective 
degree of freedom. For example, it would certainly be 
a direct interaction if the projectile bounced off a 
spheroidal nucleus, leaving it rotating. These 
rotational states are typically low lying states. 
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Compound nucleus reactions, however, involve 
many degrees of freedom due to the complicated nature 
of the nuclear excited state. It would then be 
fortunate to have a detector capable of indicating 
the number of degrees of freedom in the wave functions 
of the intercepted particles, but alas, none has ever 
been found. In a typical compound nucleus reaction 
the probability of a particle being emitted leaving 
the nucleus in a low lying state is greater than those 
leaving the nucleus in the higher excited states. 

This is due to the difficulty of barrier penetration 
which inhibits the emission of particles with energies 
just above the separation energy. Figure (1-1) shows 
a rough plot of the number of particles emitted vs. 
the excitation energy of the nucleus, where is the 
ground state energy. 




The vertical lines show the discreet spectrum. The 
detector, however, does not possess perfect resolution 
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and the energy spectrum it sees is averaged over an 
interval AS. The solid curve shows such an averaging. 
Thus it is seen that although the compound nucleus 
would rather emit particles leaving the energy of the 
nucleus near its ground state, the majority of the flux 
occurs in the higher lying states. 

Hence it is seen that inelastic scattering to 
low lying states should go predominantly by the direct 
interaction process. This theory will be developed 
in the following chapters. 



CHAPTER II 



THE SCATTERING CROSS SECTION 
Introduction 

A useful concept in relating theory to experiment 
is the differential cross section. In the usual picture 
of scattering, a collimated beam of particles impinges 
upon a target. The number of particles per second, dM, 
scattered into an element of solid angle is proportional 
to the incident flux N and the element of solid angle cl-^L 

dw = or (e>4>) W'cUi- (2-1) 

The constant of proportionality Cr( ©-,(})) has dimensions 
of area and is called the differential cross section. 

It can be a function of 0 and <j) the angular coordinates 
of dll, but is usually only a function of 0 taking the 
polar axis as the beam direction. The cross section is 
related to Schrodinger *s equation in the following manner. 
At distances far from the scattering center, the wave 
function describing the scattered particles can be 
written as the sum of a plane wave corresponding to the 
incident beam and an out going spherical wave due to the 
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scattering. 



Y - ^ ' C + 



y^-k^ef 



( 2 - 2 ) 



r 



The function 4te4>) is known as the scattering amplitude 
and contains all of the angular dependence of the 
scattered beam. 

The quantum mechanical expression for the scattered 
flux is given by"^ 




(2-3) 



where ha p is the reduced mass of the scattered system. 
Applying this to (2-2), an expression is obtained 
giving the number of particles being scattered / unit 
area / sec. 

C _ Ilk r 1 -P ( e 4*3 1 

scAr yrvp p2. 



(2-4) 



where, in the asymptotic region considered, only the 
radial component contributes significantly to the flux. 

The number of scattered particles passing through 
an element of area dA of a sphere of radius r located 
concentric with the scattering center is just dW. 
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d*4 - (s SCAr ) dA - iLkf j -P(©4 >) \ dsL. ( 2 - 5 ) 
where dn_ =r dA 

ir*^ 

The incident flux is computed using an expression 

i k-Z 

similar to (2-3), where the incident wave is £ 



S lKic = M (2-5) 

m I ■ 

The quantity clNl is then divided by the incident flux 
and the resulting expression compared with (2-1) to 
obtain the scattering cross section; 



CT(©4>) = 




(2-7) 



In the above intuitive argument, it should be 
noted that any interference between the incident and 
scattered waves has been ignored* It would thus appear 
that (2-7) is inaccurate and of limited usefulness. 

The incident beam however, has a macroscopically small 
cross sectional area and thus would only interfere with 
the scattered wave at both very small and very large 
scattering angles. Since the density of particles in 
these beams is small, they will not greatly influence 
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one another even in these cases. For a more rigorous 
treatment, the reader is directed to reference /~2 _J7, 
where it will be noted that the result is the same. 

It is the purpose of the remaining sections of 
this chapter to obtain expressions for the scattering 
amplitude, using the Schrodinger equation with an 
appropriate potential. As this is a study of the 
direct interaction, the Hamiltonian will exclude all 
flux that would go into populating the states of the 
compound nucleus and would itself describe only those 
few internal degrees of freedom pertinent to the direct 
interaction. Typically, it would describe the notions 
of the one or two nucleons to which momentum had been 
transferred. 

i 

The Schrodinger Equation 

Consider the Shrodinger equation in the center 
of mass system. The subscript f means “with respect 
to the final channel f •" 

-£,v z + * Mji 

The first two terms on the left are the kinetic energy 
and the optical interaction potential of the separating 
fragments, and is the Hamiltonian describing 

their internal motions. is the residual 



^ = E If (2-8) 
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interaction energy, with respect to the final stats 
channel, which gives rise to transitions to excited 
states of the residual nucleus. The optical potential 
is of the form 



where the term on the left is the Coulomb potential 
for a point charge 2 and an incident particle of 
charge^ . The term on the right is the phenomenological 
complex nuclear potential first introduced by Saxon 
and thus bears his name. While it will suffice here 
to give a plausibility argument for the complex nature 
of this potential, the reader is referred to and 

for a theoretical justification. The study of 
dispersion- theory in physical optics has shown that 
absorption of light in a cloudy or semi-transparent 
material can be represented mathematically by the 
imaginary term of a complex index of refraction. In 
the study of nuclear reactions, the phase of a particle 
wave impingent on a potential is shifted as it passes 
through the region of the potential and it is known 
that the scattering amplitude is completely describable 
in terms of these phase shifts.^ Completely analogous 
to the optical case, if these phase shifts are complex, 
the imag inary part would produce absorption of part of 






\r 
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the Incident beam. It is then the imaginary part 
of the complex potential which gives rise to absorption. 
This term literally absorbs out all flux which is not 
described by the potential*VO;f), such as the compound 
nucleus flux, thus eliminating it from the scattered 
beam. The denominator gives the radial dependence of 
the potential combining its short range nature with a 
diffuse edge. It will be mentioned in passing that 
many variations on this appear in the literature. 

Hov/ever, they often contain so many parameters that the 
question arises as to whether or not the resultant fits 
contain any physics or are just attempts at curve 
fitting. The variables (or. ^ in the incident channel) 
are those few internal degrees of freedom considered. 

The Hamiltonian 7 tf( is assumed to have eigen 

functions £p.) which are solutions of the equation 

H&) = hV F (^ (2-9) 

To solve equation (2-8), It is first written 
with the two potential terms isolated. 

- e) ■ if = -(u w)'* (2-10) 

To remove the ^.dependence, the equation is multiplied 
by the complex conjugate of VpC^/and integrated 
over the internal variables. 
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£ V - r 



zm 



F J 






( 2 - 11 ) 



= " <(U-vV)^d'f P 



Using (2-9) and the Hermitean nature of 
equation (2-11) can be written, using Dirac’s 
"bra-ket" notation for the integrals as 



( v 2 +£) y f = 




(2-12) 


% = 


<v F (%) \ 


(2-13) 


£ - 


2 $(e-h) 


(2-14) 



% is the coefficient in an expansion of / ^f into a 
set of eigen functions and can be thought of 

as a projection of ^if on the basis vector for the final 
channel " f The homogeneous part of Equation (2-12) 
is the Helmholtz wave equation with solutions 




± t k • r 

e — 
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General Integral Equation Solution 

Equation (2-12) is transformed into an integral 
equation through the use of a Green's function. It 
can be shown that 




is a solution to (2-12), where the Green's function 
«<£>,£) . a solution to the equation 

( 7 2+ ttp) K(tv, tf ) = - SLr~-zJ) (8-16) 

as long as , (2-16) is just the wave equation 

mentioned earlier with outgoing solution 

^ 



7 

This is expanded in a series of Legendre polynomials 

co 

=Jj* (zl * i)^ (b £) p (<^> - (2 ' 17) 

Q 

Using’ the Legendre addition theorem , (2-17) can be 
expressed in the space fixed coordinate system. 
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( 2 - 18 ) 



The primed and unprimed variables are the angular 
coordinates of and ^r^ respectively. At large 
distances from the origin, the asymptotic form 
of j^-(kr) may be used. 




( 2 - 19 ) 

( 2 - 20 ) 



We proceed by developing a series solution for the 
Green’s function. All "f" subscripts will be dropped 
until they again become necessary. The three 
dimensional Dirac delta function expressed in spherical 

Q 

polar coordinates is 

= p. £(Cos6-Cos<=>) S(W) (2-21) 



The spherical harmonics form a complete set over the 

surface of a sphere and a mathematical statement of 

0.0 

this is the following “closure relation" 
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<S(G>sS- Cos©') £C^-4>') 




( 2 - 22 ) 



The total delta function can then be expressed in' terns 
of this series. The closure relation suggests the 
following expansion for the Green* s function. 




(2-23) 



This is now inserted into the defining equation, 

(2-16), where the Laplacian in spherical coordinates is 






| 

K ^ 3\v?& 



(2-24) 



Equation (2-16) becomes 



+ ky^] 






I 



4 -i- £&» 



0 




jbo 






(2-25) 



The spherical harmonics are a solution to the following 
differential equation. 
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+ *'*')&, = O 

This is solved for the second set of square brackets 
and the txvo series are equated term by term. 

( 2 - 26 ) 

jdrilW) + (If - =---LS(f--H) (8-27) 

Lst (fj H) = ry^ 1 ) (2-28) 



The radial Green’s function is a solution to 

J. S('r-r l ) (2-29) 

Consider the homogeneous part of this equation. 

i £(krir) 4 l l< - 4M±i£) £ {kb) =0 . (2-30) 

d'r* * f* * 



As this is a second order differential equation, 
it has two independent solutions. The radial Green’s 
function, > is built from these two solutions 

by the following prescription. 
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Region I 



Region II 



G(r,r' ) 




FIGURE (2-1) 



Figure (2-1) shows a plot of some general £j(^r) 
as a function of r. It must be a solution to (2-30) 
in both regions, satisfying the boundary condition 
at the origin in region I and the boundary condition 
at infinity in region II. While it is required that 
the magnitude of ^(^/* / ) be continuous across the 
boundary at rsr', there will be a discontinuity in 
the slope. The magnitude of this discontinuity can be 
obtained by integrating equation (2-29) from r-e to 
v and taking the limit as €-^0. 









i 

r' 



(2-31) 



The second integral contributes nothing, as ^(r,r y ) 
is continuous across the discontinuity at Two 

solutions to equation (2-30) will be obtained. 
Introduce the substitution 

•{j(ki-) = ft? 



This is inserted into the homogeneous equation (2-50), 
yielding 




(2-32) 



This is Bessel’s equation of index tL + k and 




(2-33) 



where is the spherical Bessel function. 

Consideration of equation (2-17) shows this to be 
kr times the radial wave function. Here, ^(kr) is 
regular at the origin and thus satisfies the boundary 
condition in region I* The following relation exists 
between the spherical Bessel and Hankel functions.^ 

^ 4 Oc fj) 



(2-34) . 



19 



The functions and ^ Ire, respectively, the 

spherical Hankel functions of type one and type two. 

They are both solutions to equation (2-30). 

/ (') 

Asymptotically f?/ goes as an outgoing spherical wave 

j (21 

and rlg_ as an incoming one. As r approaches infinity, 

there can be only outgoing scattered waves. Hence 

a solution satisfying the boundary condition in 

, 0 ) 

region II is r}j^, kr) . The Green’s function is written 

(2 - 35) 

^ ( kr ) !->*■' 1 



subject to the condition imposed by (2-31) . 

,0) 1 if , n 






- L (2- 31a) 

, r ' 

hxK V’-V'' 

The Wronskian relation for these two functions can be 

used to evaluate the constants C]_ and Cg. 

Sturm- Liouvi lie theory’ 1 ' 2 gives the V/ronskian of two 

independent solutions to equation (2-30) as 



V^'0er),^fc»)) = A 



• (2-36) 



where A is a constant. 



20 




(2-37) 



As this is true for all r, the asymptotic forms 
of these expressions may be used in computation. 
The asymptotic form of ^(ici~) was a l rea< iy given 
in (2-19) and for large r 




Calculation of - the V/ronskian yields 

A = k 

A comparison of (2- 31a) with (2-37) shows 

c, - k^W') 

c z = fade*-') 

The Green’s function becomes 

= kr^fclrz) h^Odz) 

is the smaller of r and r' 
ty. is the greater of r and r' 



(2-39) 
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Combining equations (2-23), (2-26), (2-28), and 
(2-39), the total Green’s function is obtained. 




Equation (2- 23a) is then used in (2-15) to give 
the solution for 



The General Scattering Amplitude 



As was discussed in the introduction of this 
chapter, the scattering amplitude is the coefficient 
of the outgoing scattered spherical wave. Equation 
(2-15) will be used to find the asymptotic form 
of %. Equation (2- 23a) can be simplified for large r, 

^ ✓ /<U 

r being r, where h. can be replaced by its 
asymptotic form given in (2-39). The Green’s function 
becomes 




where the primed angular coordinates refer to K. and 
the unprimed ones to Jf p which is taken to be along 
the vector Jk p . Consider the expression in (2-18). 
Taking the complex conjugate yields 
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-Vs^r 

3 \6fJ-f _ 







Q. ~- r ~ F = 4ft, 

,^=o Jvw-£, 

whore the primed coordinates refer to Jjj.' and the 
unprimed ones to X'p . ' 

d p ^-tlc F 



K(r Fj r^ =• 



J-, 

47C/p 



e 



Using equation (2-15) 



T 

T 



c^*r i //* 

e -^e ' 



£ 









4tt^ 



ife 



r ^<Tu+wWjdW^ 



clc’V^ 

The plane wave <2 " — is a solution to the- homogeneous 

part of (2-12) and can always be added to the scattered 

solution. The scattering amplitude is then equated 

with the coefficient of -pr £ 

*> 

ft<4) = - <(v F e^ r l u+v ! 1 2 - i0) 



where the Dirac notation is again used to indicate 
integration over all internal and external variables. 

The Elastic Scattering Amplitude • 



In a description of only the elastic component 
of the scattering, the potential ^-) which 

describes transitions to excited states is set equal 
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to zero and the interpretation of the imaginary part 
of the complex potential is broadened to include 
absorption of the entire inelastic flux. The initial 
and final channels are then the same and the outgoing 
wave functionals 



y= 



(2-41) 



where V z (^ r ) and are, respectively, the 

internal and elastic scattering wave functions for 

■>/Uj 

the initial channel. satisfies the equation 




The elastic scattering amplitude is obtained from 
(2-40) 




Inelastic Scattering - Distorted V/aves 

The direct interaction inelastic scattering 
component is to be isolated from the scattered wave 
and. using techniques developed in the first part of 
this chapter, the inelastic scattering amplitude v/ill 
be found. Equation (2-8) is re-written with the 
potential which describes the inelastic transition 
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isolated. 

(-iS/+ u + ^)- E )t= -vCnr&i ^ 

Using (2-9) and the Hermitean nature of , 

the ^.dependence is removed in the same manner as 
for (2-12) 

+ = ^<Vfc|Vft^|t> (2-43) 

2^ = ^WCf» (2-44) 

/ 

Idp and ''Ip are discussed after equation (2-12) and will 
he reiterated here for the readers’ convenience. 

X — <^ P (%) | ^ > 

The homogeneous part of (2-43) has as one of its 
eigen solutions, the elastic scattering wave function 
discussed in relation to equation (2-41). A series 
solution corresponding to an outgoing v/ave will be 
obtained. Consider the homogeneous equation with 
the f subscript dropped as before. 

(y z -2i+ if P ) = O 



(2-45) 



25 



The form of the desired solution is such as to reduce 
to expression (2-18) when the optical potential is 
removed, as certainly (2-45) reduces to the Helmholtz 
wave equation. As the optical potential is a function 
of only the radial variable, ( j£,Jf) will have 

azimuthal symmetry if Jc is assumed to be along the 
z-axis* (Note expression (2-17)). The discussion 
following Equation (2-8) gave the optical potential 
to be of the form 



A plausibility argument for the complex form of the 
nuclear potential was given as well as an expression, 
typical of those appearing in the literature. No 
explicit form of Vjj is necessary in the following 
calculations. However, it will be assumed that the 
nuclear forces are short ranged and thus will 

rapidly approach zero for distances greater than the 
nuclear radius. Equation (2-45) becomes 




\p 2 + if — ~~ ^ • (2-45a) 




where 



(2-46) 
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Coulomb Wave Functions 



Consider the special case of zero nuclear 
potential. This situation would be encountered if 
the impingent beam lacked sufficient energy to 
overcome the coulomb barrier and penetrate into the 
nuclear material. The eigen solutions become coulomb 
wave functions, satisfying the equation 



The second chapter of reference /"~3_7 * s ^ excellent 
source for information regarding solutions of this 
equation. It turns out that (2-47) has a regular 
solution of the form 



If z is the polar axis, then z = r Cos© and f(r-z) has 
azimuthal symmetry. Using (2-24), the Laplacian is 



[v z +H- ZjJz] <f> c c c) = o 



(2-47) 




(2-48) 



expanded in spherical polar coordinates and the 
following substitution 's made. 




u= v-O-dos©) 
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The substitution is carried out using the chain rule. 



3-ftU) _ su 

2>x da fix 



(2-49) 



The results in the following differential equation 
for fCK) 






'•Ru)=o 



(2-50) 



The factor k is put completely within the variable 
with the transformation 



y = U u. 



(2-51) 



Equation (2-50) then becomes 




Equation (2-52) is in the form of the confluent 
hypergeometric equation which is solved in appendix A. 
A comparison of equation (A-l) and its solution 
(A-4) with (2-52) give a 

4Vi) = iFj (-Lo<|i |v) 



and 
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ee 



Fj (-Co4 j 1 1 l l<0r-2)) 



(2-55) 



A consideration of equation (A-4) shows it to be 
non-singular at the origin and therefore the solution 
<|> c is the regular coulomb wave function. The constant 
C is determined so that asymptotically 




C l < IB | 

& fkcAT 



where as usual <*> sca t is the scattered outgoing 
spherical wave. In appendix A, the following asymptotic 
solution (A- 19) is derived. 









-% £(&j& z a - J= > 

P O) 



making the substitutions, equation (2-53) becomes 




P0+ Co(j 

vcCr-i) rO-U) 



(2-54) 



the logarithm is a slowly varying function with respect 
to its argument. Hence at large distances from the 
scatterer, the first of these two terms is essentially 
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an outgoing plane wave and the second an outgoing 
spherical wave, The coefficient of the plane wave term 
was defined as being unity and thus the value of C 
is determined, 

—Tce</z . 

C = ^ rO (2-55) 

Equation (2-47) is also separable in spherical 
coordinates, thus leading to a decomposition of 
into partial waves. 

Let zb c - Egg) <§>(€►) 

. lr 

vVhen this is substituted in (2-53) two equations for 
F and © are obtained. 

5^9 +Jl(£+i) <S> = O (2-56) 

^ = (2-57) 

where ^(^-hl) is the separation constant. Equation 
(2-56) is Legendre’s equation which has Legendre 
polynomials as its solutions. 

& = ^ (Cos©) 

In equation (2-51) let 

p = - £?c kt* (2-5S) 
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With this transformation, (2-57) becomes 




Let (f = * f *(W*)**\ If) 



(2-59) 

(2-60) 



After some manipulation, the equation for V^. becomes 



+ (aA+a-/>)sMi — U+i +it<)\/ A = o 



(2-61) 



A comparison of (2-61) to (A-l) shows to be a 
solution of the confluent hypergeometric equation 



- where. 



\JjL = j FJ* l 3^3.l/0 (2-62) 

The radial coulomb wave function becomes 

Fji&f-) = A/j, e(k ffr 1 , + <■<?< I SM \-lckt-) ( 2 - 63 > 

where the Ajj_ are determined by the same asymptotic 
considerations as before, namely that behaves as an 
outgoing plane wave at large r. The asymptotic 
expression for Fjj is calculated in the same manner 
as for equation (2-53). By substitution of the 
proper parameters into (A- 19) , the following result 
is obtained. 
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A > ^ r(z*+z) 

2 X (2c) 



1 _1 

ra+i-w) 



— i (u. 1 " -&//&!& — ^0) 



P(.£-H-v <■*-). 



(2-64) 



The general theory of complex variables gives 

r(t*i ±u) = jra+i-KcOl^' Jqz - 

where Cs£ =. arg^P (j?+l (2-65) 



with this substitution, expression (2-64) simplifies to 







( 2 - 66 ) 



A consideration of the plane wave expansion 
in (2-17) shows the corresponding radial wave function 
to be j^(kr), with its asymptotic form given 
in (2-19). The constants are chosen so that these 
two functions approach each other as | r — *— o*=> 

4 _ 2^ \PCt-cUU^\€ ^ (2-67) 

A P (Zj2+ 2) ' 

The total solution for (j) c is then the sum of all the 
partial waves. This is shown by using expressions 
(2-53) and (2-55) for <j> c and the contour integral 




representation of the confluent hyper geometric - 
function given in (A-7) 




* [<a Ur ' 



Q-e^n-y 



(2 - s8> 



using equation (2-17), the exponential e ikz ^ 1 ** t ^is 
expanded into a series of Legendre polynomials 
where z = r CosO 

<j> c - y l&fajL+i) ^Jjc.h) f£ (Cot.®) (2-69) 

jZrO 



V/here 












(2-70) 



a consideration of equations (2-30), (2-33), (2-57), 
(2-63), and (2-67) with equal to zero in the last 
three, gives 



ff(&4 1-^tx) (2-7i) 



using the series form of the confluent hypergeometric 
function (A-4), equation (2-71) becomes 




-? S P LK X X V ~ P(i*i+s)(-aix ) 
sTo r (2^+24 s) si 



(2-72) 



Letting x =: kr(l-t), (2-72) is substituted into 
(2-70) and the order of summation and integration 
is reversed. This gives 



hCuo = 0*1* 

a o -'&***) pc-w). 






x rfefrHsj 

^7 d R(2&*2+S) 





(2-73) 



The integral given in the parenthesis is a contour 
integral representation of the beta function as given 
in expression (A-6) 




-O'-** 0 *) ' P(^)p(o-s^u ) 
pU+s+1) 



with this substitution, (2-73) becomes 






y y~ p 6<?+a) (oc^) 

* ^To ^Ca^a-fs) p(£-w4-c©<) Sj 



By equation (A-4), this is written 



Of 



X Te^ fC * +U ^ +al ~^ 



(2-74) 



(2-74) is compared with (2-63) and (2-67), giving 



=■ , 5 > ^ 

■ ki~ 






(2-75) 



where is defined in egression (2-65). Expressions 
(2-69) and (2-75) are combined and the partial wave 
expansion of <j> c is obtained 



j> c = - 1 - y^(*M)e^ % S) (2-76) 



Optical Wave Functions 

The eigen solution to equation (2-45a) corresponding 
to an outgoing wave will now be obtained. It was shown 
in the previous section on coulomb waves, that if the 
nuclear potential is neglected, the equation is 
separable in spherical polar coordinates, and since Vjj 
is a function of only the r variable, the separability 
is unaffected. If the Laplacian is expanded and the 
variables separated, the angular wave function is 
unchanged. The equation for the radial coordinate is 




2 



— = 0 (2-77) 



The total solution is 



/^=> 



X m ( VW -p2. Qt <%(&■)%&**) t2_78) 



where, as usual, the constants are determined 
from the asymptotic considerations. 

It \va3 3ho\vn in equation (2-64) that the regular 
asymptotic form of the radial coulomb wave function 
could be built from a linear combination of the two 
asymptotic exponential solutions 



where the plus and minus signs correspond to outgoing 
and ingoing waves, respectively* As the nuclear 
potential was postulated as being short ranged, these 
may also be used in determining the asymptotic form 
of ^(kr) . This is seen by substituting the 
exponentials into equation (2-77) with the following 






aQnZkV + <%) 

e x 



result 
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. __ ±t(jcr-edUakjt~- 

C^+oi^e 




O 



It is seen that as long as Vjt approaches zero at least 
as fast as r~ 2 , an asymptotic • form of j^,(kr) can also 
be formed from these exponentials with error certainly 
no greater than that inherent in the coulomb solutions. 
The phenomenological nuclear potential postulated in the 
earlier, part of this chapter meets this requirement. 

From consideration of the corresponding coulomb solution, 
the regular asymptotic solution to the radial optical 
wave function is 

SinCkir- + (2-79) 

where <$^is an additional phase shift which the wave 
would experience in passing through the region of the 
nuclear potential* The constants Cj^ are now determined 
from the following two considerations: (1) The nuclear 

potential, being short ranged cannot affect the incoming 
asymptotic solution and, (2) VJhen % is set equal to 
zero must reduce to the coulomb wave. As it was 
shown that the solution could be written as a linear 
combination of the ingoing and outgoing coulomb waves , 
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the first requirement becomes 

Ajt, Sin Ck; It - —c^-2n k.\r — O-Tp -k<T^ -fr 

"S. 

-ICi ct — ^j2n2ich — Jtcr-t- cn) 

^ 21 (2- so) 

0 L(jcl~-*o£<2^2/cr-J 2ST-t<^:') 

"+• Djz. HEi *• 



The solution of this equation yields 




(2-81) 

(2-82) 



In order for (2-78) to reduce to equation (2-76) 
when V N and ihus are zero. 



cj =z jL? && +i) Q °*- 

he 



(2-83) 



and the total wave function becomes 

y^=- £(<«&) (2-84) 

G > *]'32. “* (To, Hr 



\7here 



(2-85) J 



Using the Legendre addition theorem to ivrite this 
in terms or the angular coordinates of the vectors 
k and r, equation (2-84) becomes 



The following development parallels that used 
in the study of the general integral equation, the main 
difference being that distorted waves are used instead 
of plane waves. The inelastic component of the 
scattering has already been isolated in equation (2-43) 
and the corresponding integral equation for is 



following equation. Note that the f subscripts have 
been dropped as before. 




(2-S6) 



where the primed variables will be taken as the 
coordinate corresponding to k. 



Inelastic Integral Equation 




^■^(r^ 1 ), the Green’s function is a solution to. the 



(v 2 -^ ii)J?i '&:) - - SCjc-ji) 



( 2 - 88 ) 



It was seen in equations (2-21) and (2-22) 
that the Dirac delta function could be expressed 
as a series of spherical harmonics as follows 



Equation (2-88) differs from (2-16) in the previous 
treatment only in the addition of a potential term 
which contains no angular dependence. The series 
expansion of the Green's function defined by (2-88) 
will then have the same angular dependence as (2-23). 
and can be written 



When this is inserted into (2-88) using the Laplacian 
given in (2-24), the radial Green's function is seen 
to be a solution of the equation 



fC h8Q -MM ) (2 _ 90) 



The reader is referred to equation (2-25) and its 
subsequent discussion, where more detail is given 
to a similar calculation. 







/-a 




Letting 




(2-91) 
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equation (2-90) becomes 

(sjl ' + ic -v. (r) ') (2-92) 

The radial Green’s function is obtained from two 
independent solutions to the homogeneous part of this 
equation, one regular at the origin, and the other 
corresponding to a purely outgoing scattered wave 
satisfying the boundary condition at infinity* The 
regular solution ^^iis obtained from equation (2-77) 
with its asymptotic form given in (2-79). 

Let3£(J'/}be the solution to ( 2 - 777 ) and thus to the 
homogeneous part of (2-92) with outgoing boundary 
conditions. It was shown in equations (2-80) and (2-81) 
that the asymptotic regular solutions could be written 
in terms of the two irregular exponential solutions 

drc(^ 11 "- 

e 

where the plus and minus signs correspond, respectively, 
to outgoing and incoming waves. 

Thus, asymptotically, J3^(k,r) becomes 

±C (kt — Art + #}%) 

r * 



(2-93) 



The Green's function is written 



gw 



r 

£ z H x (k>-) 



lr< tr* 

r? h* 



(2-04) 



where, analogous to equation (2-31a), the discontinuity- 
in the derivative of ^(r^r') evaluated at r' is written 

ezifa (kh) ~ (2-ss) 

hah' hr 1 

The constants G-^ and are determined from the 
V/ronskian relation 



WC^ j^jl) - -A' (2-96) 

where A* is a constant. The Wronskian is evaluated 
using the asymptotic forms of and and f ound 
to be 



XI d U 





(2-97) 



A comparison of equation (2-97) with (2-95) yields 
the values of the constants. 



rv Mu ilcr > 




(2-98) 
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Using the notation of (2-39), the radial Green's 
function becomes 






~ (2- 94a) 



and the total Green’s function is obtained from 




(k%) 



(2-89a) 



This is then used in expression (2-87) to obtain the 
outgoing scattered wave. 



Inelastic Scattering Amplitude 

Using (2-87), an expression f or ^in the asymptotic 
region is sought* For large r, r >>r » , Jfy ( kh ) can be 
replaced by its asymptotic form given In (2-93) and 
the Green’s function written. 

(2-99) 

X ** 

where OV^is defined in expression (2-85) and r^ is • 
taken along the direction of k^.- . •, 

Consider the following part of (2-99) 
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( 2 - 100 ) 



Upon comparison with (2-86), those are recognized 
as being similar. 

Expression (2-100) is identified with the complex 
conjugate of the ;"time reversed" optical wave function, 
and is denoted by^C.(kr). In the "time reversed" 
scattering problem, the initial and final channels 

y * 

have exchanged roles andX-Ykr-) is a solution to 
equation (2-45a) with incoming boundary conditions, 
where has been replaced by its complex conjugate 
to simulate flux re-inserted into the beam. ‘The 
function is perhaps best generated using the previous 
calculations for^^(k,r) and the Wigner Relation^ 






( 2 - 101 ) 



Using equation (2-86), 

yrc&) = 

14 

The parity of the spherical harmonic is 




( 2 - 102 ) 
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Hence, equation (2-102) becomes 







(2-103) 



'.which agrees with (2-100) 

The Green’s function, (2-99), is then re-written 
incorporating. the "time reversed" wave function. 



Mi 







* x%$> 



(2-99a) 



Using equations (2-87) and (2-99a), the asymptotic 
egression for "%■ is obtained. 







_e 












The first term is the elastic scattering or optical 
wave function which is a solution to the homogeneous 
equation (2-45) and is necessary should the initial 
and final channels coincide. This would be the case 
of elastic scattering. The inelastic scattering 
amplitude is the coefficient of the outgoing spherical 
wave . 



I V I > (2-105) 



With the inelastic component of the scattering 
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amplitude separated, the total scattering amplitude 
can be written, not only as in equation (2-40), 
but also as a sum or its elastic and inelastic 
components. Using equations (2-42) and (2-105), 
the total scattering amplitude is written. 



t areif 






(2-106) 



Matrix Element Formalism 



Customarily, the scattering amplitude is not 
used in the literature, but rather the matrix element 
giving the probability of transition from an initial 
to some final scattering state. A consideration of the 
expressions for the scattering amplitude shows the 
matrix element to be given by 

r FI = - (2 - io7 > 

where the subscripts i and f refer to the initial 

and final scattering states, respectively. 

In this formalism, equation (2-106) becomes 

15 

the Gell-Mann- Go ldberger relation . 

u +<Xf-^lVl (2-108) 

which, while derived here for the specific cs.se 
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of the optical wave functions, is quite general and is 
usually obtained from a unitary transformation of the 
general scattering amplitude matrix element. 



» / 

t fx = <y F e ~ F ~' r t u*v| 



(2-109) 



where u. and are any tv/o potentials that have been 
isolated. 

The angular distribution or cross section is 
written in this formalism using equation (2-7) 



OTei) = 



iZTTk*) 7 - 




( 2 - 110 ) 



The summation sign indicates a twofold operation; 

(1) due to the inability of the detector to distinguish 
betv/een spin orientations, the projections of the 
final state spins are summed, and (2) the projections 
of the initial state spins are averaged as when time 
and thus the direction of particle motion is reversed, 
the direction of any spin must also be reversed to 

1 g 

preserve time reversal invariance. It is seen that 
if these projections are averaged, producing "an -"average 
spin" initial state wave function, the problem is 
averted. A suitable modification would have to be 
made, however, if polarization experiments are to be 



considered 
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It now only remains to specify a suitable potential 
and scattering wave function and, using (2-110), the 
cross • section can be obtained. This problem is 
considered in the following chapter. 
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CHAPTER III 



EVALUATION OP THE DIFFERENTIAL CROSS-SECTION' 

FOR COLLECTIVE EXCITATIONS 
Introduction 

An exact calculation of the matrix element T^ 
is at best, a difficult problem. Transforming 
Schrodinger * s equation to its integral equation 
counterpart has not given its solution, but has only 
provided a convenient way to obtain the scattering 
cross section once the wave function has been found. 

The specific interaction potential is, in general, 
not known and thus an exact solution for ^ is not 
easily found. For example, if the angular momentum 
transferred in the reaction is non- zero, the interaction 
potential will not, in general, be a scalar function 
of position and normal methods for solution of the 
Schrodinger equation fail. Physically meaningful 
results, however, may be obtained with appropriate 
approximations for the wave functions and the 
interaction potential. It is the purpose of this 
chapter to discuss the techniques involved, and to 
show their applicability by fitting experimental 
angular distributions* 

In each of the expressions for T^, four 
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integration variables are employed. They are 
j*i, ^ and S~ £ where each is a vector quantity. 

Undoubtedly, both the initial and final position 
variables, and the initial and final internal variables 
are related, but the exact nature of this relation 
is unknown. It would be convenient to have but one 
position variable and one internal variable when 
performing these integrations, hence for want of this 
exact relation let 




= Sf = 5 



This is called the zero-range approximation. 

The angular* distribution is now evaluated using 
the general scattering matrix element. 

Tpx s U+v I (3-2) 



It will be assumed that the wave function can be 
at least partially factored and takes on the form 






( 3 - 3 ) 



where the wave function is the eigen function 

for the initial nuclear state, and ^^.r > ^ ) is a 
solution to the equation 
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(3-.4) 



It should be noted that the Hamiltonian describing 



variables can be treated as constants and the function 
/'j; depends on them only in a parametric sense. 

is then obtained by a two step calculation; 

(1) equation (3-3), a function of only the r variable, 
is solved and the resulting solution used to compute 
the transition amplitude for a fixed value of , 



and (2) this result is multiplied by v^(^) and the 
projection of their product on the final state f is 
f ound. 

Physically, the above corresponds to the assumption 
that the initial nuclear state described by the function 
v i (i^) changes so slowly with respect to the transit 
time of the projectile, that the scattering wave function 
is completely determined by the value of Vj^ir?) at the 
instant of the reaction. This is analogous to “stop 
motion" photography, where a moving object has its 
motion frozen on film through the use of a short duration 



the internal motions is not present. Hence the 



where r Fx (§) =• u+i/| 9^ } > 



(3-5) 



53 



flash of light. In the time considered, the 
variables are almost stationary and thus it is seen 
that tho Hamiltonian, , can be neglected in. the 

integral equation leading to (3-2) and in (5-3). 

These assumptions are most valid when high energy 
projectiles are used and when the nuclear excited 
states are collective states. Since the excitation 
energy must be shared with the fairly large mass 
of the nucleus, the collective variables exhibit a 
slow change. 

The adiabatic theory (Diffraction Scattering) 
will now be developed and, using the computed angular 
distributions, will be compared with experimental 
alpha particle scattering data. 

Diffraction Scattering 
The Fraunhofer Approximation 

The calculation is most conveniently made in the 
scattering amplitude formalism; hence (3-5) becomes 




Equation (3-3) is multiplied through by 2m/h 2 and the 
resulting expression is substituted into (3-6) thus 
transforming away the dependence on the interaction 
potential. 



